Modelization of Thermal Fluctuations in G Protein-Coupled 

Receptors 

C. Pennetta 1 , V. Akimov 1 , E.Alfinito 1 , L. Reggiani^G. Gomila 2 , G.Ferrari 3 , L. Fumagalli 3 and 

M. Sampietro 3 

1 Dipartimento di Ingegneria dell'Innovazione, Universitd di Lecce and 
National Nanotechnology Laboratory-INFM, Via Arnesano, 73100 Lecce, Italy. 

2 Department d'Electronica and Research Centre for Bioelectronics and Nanobio science, 
Universitat de Barcelona, C/ Josep Samitier 1-5, 08028 Barcelona, Spain. 

3 Dipartimento di Elettronica ed Informazione, Politecnico di Milano, P.zza Leonardo da Vinci 

32, 20133 Milano, Italy 

Abstract 

We simulate the electrical properties of a device realized by a G protein coupled re- 
ceptor (GPCR), embedded in its membrane and in contact with two metallic electrodes 
through which an external voltage is applied. To this purpose, recently, we have proposed 
a model based on a coarse graining description, which describes the protein as a network 
of elementary impedances. The network is built from the knowledge of the positions of the 
C a atoms of the amino acids, which represent the nodes of the network. Since the elemen- 
tary impedances are taken depending of the inter-nodes distance, the conformational change 
of the receptor induced by the capture of the ligand results in a variation of the network 
impedance. On the other hand, the fluctuations of the atomic positions due to thermal 
motion imply an impedance noise, whose level is crucial to the purpose of an electrical de- 
tection of the ligand capture by the GPCR. Here, in particular, we address this issue by 
presenting a computational study of the impedance noise due to thermal fluctuations of the 
atomic positions within a rhodopsin molecule. In our model, the C Q atoms are treated as 
independent, isotropic, harmonic oscillators, with amplitude depending on the temperature 
and on the position within the protein (a- helix or loop). The relative fluctuation of the 
impedance is then calculated for different temperatures. 

1 Model and Results 

G protein-coupled receptors (GPCRs) constitute the largest family of trans-membrane receptors, 
with functions going from revealing light and smells to the individuation of drug and virus in- 
truders [1]. For this reason, many efforts are devoted to the study of their properties [1]. In 
particular, we are interested to develop electronic nanobiosensors based on GPCRs. Actually, the 
detection of an electrical signal from hybrid nanodevices based on a single or few receptors and 
associated with the capture of the ligands, is a challenging goal, rich of potential applications 
[2, 3]. Here, our aim concerns with the calculation of the electrical properties of a device realized 
by a G protein coupled receptor [1], embedded in its membrane and in contact with two metallic 
electrodes through which an external AC voltage is applied [4]. To this purpose, recently, we 
have proposed a model based on a coarse graining [5] description, which describes the protein as 
a network of elementary impedances [4]. The network is built from the knowledge of the position 
of the C a atoms of the amino acids [6], which are taken as the nodes of the network [7]. At 
least in the case of rhodopsin (photonic receptor) these positions are known by X-rays diffrac- 
tion experiments [6] for both, the basic and the most stable excited state (metarhodopsin) [6]. 
Though these positions are generally unknown for the other receptors, a coarse graining, complex 
network approach offers the possibility of taking advantage of the common topology of the GPCR 



family [1]. In fact, all GPCRs share a seven-helices trans-membrane structure, where the seven 
a-helices are interconnected by extracellular and intracellular loops [1, 6]. Additionally, there are 
two terminal chains: an extracellular chain (N terminus) and an intracellular chain (C terminus) 
[1, 6]. We assume that the amino acids interact electrically among them and that charge transfers 
between neighboring residues and/or changes of their electronic polarization [8] affect these inter- 
actions [4]. Accordingly, a link is drawn between any pair of nodes neighboring in space within a 
given a distance, d = 2R a (R a electrical interaction radius) [5, 7] and an elementary impedance 
is associated with each link [4]. Moreover, two extra nodes can be introduced in the network, 
associated with the electrodes, which are linked to a given set of amino acids, depending on the 
particular geometry of the contacts in the real device. The elementary impedance is taken as the 
impedance of a RC parallel circuit (the most usual equivalent passive AC circuit) [4]. Precisely, 
by denoting with Z it j the impedance associated with the link between the i-th and j-th nodes, 
separated by a distance kj, we take [4]: 

where u> is the frequency of the external voltage, p the resistivity of the resistor, e the vacuum 
permittivity and e it j the relative dielectric constant of the capacitor [4] expressed in terms of the 

intrinsic polarizabilities ai and Oij of the corresponding amino acids [8]. By taking the values 

R a = 12.5 A and p = 10 9 Qm, we have found [4] that the conformational change of the receptor 

induced by the capture of the ligand (i.e. the transition rhodopsin — > metarhodopsin) implies 

a significant variation of the impedance at all frequencies, and in particular we have found a 

variation of about 20 % in the static value of Re[Z] [4]. On the other hand, the fluctuations 

of the atomic positions due to the thermal motion [5, 9, 10] imply an impedance noise, whose 

level, in comparison with the impedance change due to variation of conformation and with the 

electrode/amplifier noise, is crucial to the purpose of an electrical detection of the ligand capture 

by the GPCR. Therefore, here we consider the effect of the thermal atomic motion on the electrical 

response to an external field of a rhodopsin molecule. To this purpose, we allow the nodes of the 

network (C a atoms) to fluctuate around their equilibrium positions. For the sake of simplicity 

and to get a qualitative estimation, we describe the system of coupled oscillators as a set of 

independent, isotropic, harmonic oscillators. When the oscillators are in their ground state, their 

positions, r, referred to the equilibrium ones, are distributed with a probability density: 



where M is the average mass of the amino acids, u = y*y/M the oscillator frequency, 7 the 
force constant and < x 2 >= 1/3 < r 2 > is the mean square displacement of the oscillator from 
its equilibrium position along the x-direction. If each oscillator is in contact with a thermal bath 
at temperature T, the value of < x 2 > at the equilibrium is: 



1 k B 9 k B 9 



< X >_ 2 7 + 3 7 exp[9/T] - 1 (3) 



where 9 = huj /k B . When T 3> 9, Eq. (3) simplifies in: 



< x 2 >tt k B T/3-f (4) 
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Figure 1: Simulation of the modulus of the network impedance, \Z\, versus time at T = 300 
K. The time is expressed in simulation steps while \Z\ has been normalized to its average value 
(shown by the gray line). 



Figure 2: Relative root mean square fluctuation of the impedance modulus as a function of the 
temperature (in K). The three set of data refer to simulations performed wit the ratio of the helix 
and loop force constants, 71/72, equal, respectively, to 1 (full circles), 50 (open squares) and 60 
(full diamonds). The dashed curves show the best-fit with an exponential function. 
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Of course, for an arbitrary temperature T ^ the wave function of the oscillator is a super- 
position of several excited states and the density probability cannot be expressed in the simple 
form of Eq. (2). However, again for simplicity, we keep this expression and we account for the 
effect of the temperature by assuming in Eq. (2) < x 2 > given by Eq. (4). Moreover, as a first 
approximation, we take the mass M and the force constant 7 equal for all the oscillators [10], 
and precisely: M = M = 100 Dalton, 7 = 2.5 KJ mole" 1 A" 2 [5]. This choice provides 9 = 12 
K. Thus, the condition T 3> 9 is satisfied at room temperature. Figure 1 shows the results of 
simulations at 300 K of the modulus of the network impedance, \Z\, versus time (the modulus has 
been normalized to its average value). On the other hand, it is well known [1, 9, 11] that loops 
and terminals are very flexible structures compared with the quite rigid a-helices. Therefore, to 
overcome the crude approximation of a unique force constant for all the oscillators, we consider 
two different spring elastic constants, 71 and 72, for oscillators belonging to the a-helices and to 
loops/terminals, respectively, with 71 > 72. Figure 2 shows the relative root mean square fluctu- 
ation of the modulus of the impedance, i.e. the root mean square fluctuation, a, normalized to 
the average value of \Z\, as a function of the bath temperature. The three sets of data report the 
results of simulations performed by taking the ratio, F = 71/72, equal 1 (full circles), 50 (open 
squares) and 60 (full diamonds). In all cases the values of 71 and 72 have been chosen to keep 
constant the average value of the force constant, 7 — (V-W)Z)7i) where the sum is performed 
over the whole number N of considered oscillators. The dashed curves in this figure represent 
the best-fit with exponential functions. We can see that once the higher flexibility of loops and 
terminals is accounted for, the relative fluctuation of the impedance increases significantly and it 
becomes strongly sensitive to the temperature. ^From this study we can conclude that a careful 
modelization of thermal motion [9, 11] is necessary to provide reliable estimates of the relative 
fluctuation of the impedance. 
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